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0. INTRODUCTION 

The so-called nonlinear stationary-phase-steepest-descent method for the as- 
ymptotic analysis of Riemann-Hilbert factorization problems has been very suc- 
cessful in providing 

(i) rigorous results on long time, long range and semiclassical asymptotics for 
solutions of completely integrable equations and correlation functions of exactly 
solvable models, 

(ii) asymptotics for orthogonal polynomials of large degree, 

(iii) the eigenvalue distribution of random matrices of large dimension (and 
related universality results), 

(iv) proofs of important results in combinatorial probability (e.g. the limiting 
distribution of the length of longest increasing subsequence of a permutation, under 
uniform distribution) . 

Even though the stationary phase idea was first applied to a Riemann-Hilbert 
problem and a nonlinear integrable equation by Its ([I], 1982) the method became 
systematic and rigorous in the work of Deift and Zhou ([DZ], 1993). As a recognition 
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2 SPYRIDON KAMVISSIS 

of the fruitfulness of the method, Percy Deift was invited to give a plenary address 
to the recent ICM in Madrid, on the subject of universality in mathematics and 
physics. Of course, the main mathematical tool in proving universality theorems 
has been the nonlinear stationary-phase-steepest-descent method. 

In analogy to the linear stationary-phase and steepest-descent methods, where 
one asymptotically reduces the given exponential integral to an exactly solvable 
one, in the nonlinear case one asymptotically reduces the given Riemann-Hilbert 
problem to an exactly solvable one. 

Our aim here is to clarify the distinction between the stationary-phase idea 
and the steepest-descent idea, stressing the importance of actual steepest-descent 
contours in some problems. We claim that the distinction partly mirrors the self- 
adjoint / non-self-adjoint dichotomy of the underlying Lax operator. To this aim we 
first have to review some of the main groundbreaking ideas (due to Percy Deift and 
his collaborators) appearing in the self-adjoint case; then we describe recent results 
([KMM], [KR]) in the non-self-adjoint case, that we see as a natural extension. We 
mostly use the defocusing / focusing nonlinear Schrodinger equation as our working 
model, but we also digress to the KdV at some point. 

We stress both here and in the main text that an extra feature appearing 
only in the nonlinear asymptotic theory is the Lax-Levermore variational problem, 
discovered in 1979, before the work of Its, Deift and Zhou, but closely related to 
the so-called "g- function" which is catalytic in the process of deforming Riemann- 
Hilbert factorization problems to exactly solvable ones. 

1. THE LINEAR METHOD 

Suppose one considers the Cauchy problem for the linearized KdV equation: 
Ut + u xxx = 0. It can of course be solved via Fourier transforms. The end result 
of the Fourier method is an exponential integral. To understand the long time as- 
ymptotic behavior of the integral one needs to apply the stationary-phase method 
(see c.g.[E]). The underlying principle, going back to Stokes and Kelvin, is that the 
dominating contribution comes from the vicinity of the stationary phase points. 
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Through a local change of variables at each stationary phase point and using inte- 
gration by parts we can calculate each contributing integral asymptotically to all 
orders with exponential error. It is essential here that the phase 1 x£ — £ 3 t is real 
and that the stationary phase points are real. 

On the other hand, suppose we have something like the Airy exponential inte- 
gral 

1 f°° s 3 

(1) Ai(z) = — / cos(— + zs)ds, 

n Jo 3 

and we are interested in z — > oo. Set s = z 1 / 2 t and x = z 3 / 2 . 

r l/3 roo ,3 

(2) Ai(x 2 / 3 ) = — J exp(tx(- + t))dt. 

The phase is h(t) = ^ + t and the zeros of h'(t) = (t 2 + 1) are ±i. As they arc not 
real, and since the integrand is analytic, one must deform the integral off the real 
line, along particular paths. These are referred to as steepest descent paths. They 
are given by the simple characterization 

(3) Imh(t) — constant. 

In our particular example, the curves of steepest descent are the imaginary axis 
and the two branches of a hyperbola. By deforming to one of these branches, we 
finally end up with with integrals which can be analyzed directly, using the so-called 
Laplace's method (which is simpler than the stationary phase method). We thus 
recover asymptotics valid to all orders. 

The nonlinear method generalizes the ideas above, but also employs new ones. 

2. THE NONLINEAR METHOD 

(i) The Stationary Phase Idea 

Consider the defocusing nonlinear Schrodinger equation 

idtip + d\i) - |VfV = o, 

(4) 

under ip(x,0) = tpo(x), 

J £ is the spectral variable 



4 SPYRIDON KAMVISSIS 

where the initial data function lies in, say, Schwartz space. The analog of the 
Fourier transform is the scattering coefficient r(£) for the Dirac operator 

L _ ( id x iipo{x) \ 
\ -iipo( x ) ~ id x ) ' 

Suppose we are now interested in the long time behavior of the solution to (4). The 

inverse scattering problem can be posed in terms of a Riemann-Hilbert factorization 

problem. 

THEOREM. There exists a 2x2 matrix Q with analytic entries in the upper 
and lower open half-planes, such that the normal limits Q+,Q-, as £ approaches 
the real line from above or below respectively, exist and satisfy 

and Hm^ooQiO = I- 

The solution to (4) is recovered via 

(6) iP(x,t) = -2Km^oo£Qi2(0- 

It was first realized by Its [I, IN], that the leading order behavior of the long time 
asymptotics for the solution of (4) can be described by replacing the problem (5) 
by a "local" model Riemann-Hilbert problem located in a small neighborhood of 
the stationary phase point £ = — satisfying 0'(£o) = where 9 = £x + 2£ 2 t, 
or, equivalently to a problem similar to (5) but where £ is replaced by the constant 
value £o (see (7) below). But no idea was given on how to show that this solution 
of the full problem and of the model problem are actually close to each other. To 
show how to do this, was the work of Deift and Zhou. 2 The basic ideas of [DZ] 
have been used in all works on the stationary-phase-steepest-descent-method since. 
They include: 

1. Appropriate lower/diagonal/upper factorizations of jump matrices. 

2. Equivalence of the solvability of inhomogeneous matrix Riemann-Hilbert 
problems to the invertibility of associated singular integral operators. This idea 



2 See [DZ] for a comprehensive review of the history of the problem before the use of Riemann- 
Hilbert problem tecniques. 
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goes at least back to Gohberg [CG]. A crucial contribution of Bcals and Coifman 
[BC] was to make this precise for contours with sef-intersections (which they needed 
in their study of the inverse scattering of first order systems). These ideas where 
further developed by Zhou [Z] who provided a very useful existence theorem for 
matrix Riemann-Hilbert problems with jumps and jump contours satisfying some 
special Schwarz reflection type symmetries and an integral formula expressing the 
solution of the Riemann-Hilbert problem in terms of the inverse of a particular 
weighted Cauchy operator depending on a given factorization of the jump matrix 
and thus taking advantage of the factorization mentioned above. Perturbing Zhou's 
formula provides a nice way to show that under some conditions, small changes in 
the jump data result in small changes in the solution. 

3. Introduction and solution of auxiliary scalar problems. 

Following analyticity and the above ideas one ends up with a problem on a 
small cross centered at the stationary phase point. Using a rescaling the Riemann- 
Hilbert problem is rcscalcd to a new problem on an infinite cross. After deforming 
the components of the cross back to the real line, it is equivalent to the following 
problem on the real line: 

l-K£o)| 2 -r*(i )\ _,.. t2 



where v is a constant depending only on £ and cr 3 = ^ ^ ^ j is a Pauli matrix. 

( l-k(f)| 2 -r*(£)\ 
So the jump matrix I (£) 1 ) °^ or ^> ma l problem is replaced 

by its value at £o- 

Problem (7) can be solved explicitly. Written in terms of the new unknown 
H{^)exp{—i^ 2 a 3 ) it has a constant jump and can thus be reduced to a first order 
linear matrix ODE ([I]). 



(ii). The finite-gap g- function mechanism and a "shock" phenomenon with no 
linear analogue. 
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An important step halfway between the leap from the "stationary phase" idea 
to the general definition of a " steepest descent contour" is the introduction of the 
so-called finite-gap g-function mechanism. The g-function was introduced in [DZ95] 
in the special case of genus and in [DVZ94] in the special case of genus 1 but the 
full force of the finite-gap g-function idea and the connection to the Lax-Levermore 
variational problem was first explored in the analysis of the KdV equation [DVZ97] 

u t - 6uu x + e 2 u xxx = 0, 

u(x,Q) = u (x), 

in the limit as e — > 0. Assume for simplicity, that the initial data are real analytic, 
positive and consist of a "hump" of unit height. 
The associated RH problem is 



S+(z) = S-(z)[ L; + i„3 t ),lmz = 0, 

\ r{z)e « 1 J 

and lim z ^ 00 S(z) = (1,1), 
where r is the reflection coefficient for the Schrodinger operator with potential uq. 
The solution of KdV is recovered via 

d 

u(x,t;e) = —2ie—S 1 (x,t;e), 

where Si is the resdiue of the first entry of S at infinity. The reflection coefficient 
r also depends on e. In fact, the WKB approximation is 

r(z) ie - X[o,i]{z) 

l-\r(z)\ 2 ~ e • , 

where 

p(z) = x + z+ / [z - (z 2 - u (x)) 1/2 ]dx, 

J x + 

t(z) = Ref iMx) - z 2 ) 1/2 dx 
and x+(z) is the largest solution of uq(x + ) = z 2 . 

[DVZ97] introduce the following change of variables S(z) — S(z)e 9< ^ 3 where 
g is a scalar function defined by the following conditions. 

1. g is analytic off the interval [0, 1], the normal limits g+,g~ of g exist along 
[0, 1] and g vanishes at infinity. 
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2. "Finite gap ansatz". Define h(z) = g+(z) + g~(z) - 2p + Atz 3 + xz. There 
exists a finite set of disjoint open real intervals ("bands") Ij £ [0, 1] such that 

3a. For z e Ujlj, we have — t < (g + — g~)/2i < and h' = 0. 

3b. For z 6 [0, 1] \ Ujlj, we have 2ir = g + - g_ and h' < 0. 

The conditions above are meant to determine not only g but also the band- 
gap structure in [0, 1]. In general (for any data uq) it is not true that the above 
conditions can be satisfied. It is believed however that under the condition of 
analyticity a g- function satisfying the "finite gap ansatz" exists. (In fact [K00] 
gives a proof of the "finite gap ansatz" in the analogous problem of the continuum 
Toda equations.) Assuming that there is a g- function satisfying the three conditions 
above one can show that the RH problem reduces to one supported on the bands 
Ij with jumps of the form 



and in fact, because of (iib), h(z) is a real constant on each band Ij. This RH 
problem can be solved explicitly via theta functions. The details in [DVZ97] involve 
the so-called "lens" -argument: auxiliary contours are introduced near pieces of the 
real line (one below and one above each band/gap) and appropriate factorizations 
and analytic extensions are used, very similarly to the subsection above. The 
conditions for g above are chosen precisely to make the lens argument work. 

As is remarked in [DVZ94] the fact that the new RH problem is on slits "is a 
new and essentially nonlinear feature of our nonlinear stationary phase method" . 
Unlike [DVZ94] where it is defined explicitly via an integral formula, in [DVZ97] 
the g-function is only defined implicitly via the conditions above, which may or 
may not admit a solution. 

(iii) The Lax-Levermore Variational Problem 




The g-function satisfying conditions (i), (ii), (ha), (iib) can be written as 
g{z) = / log(z - r])d(i(ri) 
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where u is a continuous measure supported in Ujlj- In a sense, the reduction of 
the given RH problem to an explicilty solvable one depends on the existence of a 
particular measure. Conditions (i), (ii), (iia), (iib) turn out to be equivalent to a 
maximization problem for logarithmic potentials under a particular external field 
depending on x,t, uq(x) over positive measures with an upper constraint. This is 
related to the famous Lax-Levermore Variational Problem [LL]. Even though it 
appears as an afterthought in [DVZ97] (though clearly serving as inspiration), it 
seems that its analysis is essential for the justification of the method (as in [K00]). 

3. STEEPEST DESCENT CONTOURS 

Having reviewed some essential ideas in the previous sections, we are ready to 
consider the focusing NLS equation, following [KMM]. 

ihdt*l> + y^V + IVlV = o, 

under data ip(x,0) = iio(x)- 

Note that the Lax operator 

ihd x -iip a (x)\ 
-iipo(x) ~ihd x J ' 

is non-self-adjoint. We shall see that the deformation of the scmiclassical RH prob- 
lem can be no more confined to a small neighborhood of the real axis but is instead 
fully two-dimensional. A steepest descent contour needs to be discovered! 3 

For simplicity consider the very specific data ipo (x) = Asechx where A > 0. 
Let x-(n) < x + (n) be the two solutions of sech 2 (x) + rj 2 = 0. Also assume that 
h = A/N and consider the limit N — > oo. It is known that the reflection coefficient is 
identically zero and that the eigenvalues of L lie uniformly placed on the imaginary 
segment [— iA, iA}. In fact the eigenvalues are the points Xj = ih(j + 1/2), j = 
0, N — 1 and their conjugates. The norming constants oscillate between —1 and 
I. 

3 By the way, in the long time asymptotics for the above with 6 = 1 a collisionlcss shock 
phenomenon is also present; for x, t in the shock region the deformed RH problem is supported 
on a vertical imaginary slit. (See [K96].) But here, we rather focus on the semiclassical problem 
h — > which is far more complicated. 
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The associated RH problem is a meromorphic problem with no jump: to find a 
rational function with prescribed residues at the poles Xj and their conjugates. It 
can be turned into a holomorphic problem by constructing two loops, one denoted 
by C say, encircling the Xj and one, C* , encircling their conjugates. We redefine 
the unknown 2x2 matrix inside the loops so that the poles vanish (there is actually 
a discrete infinity of choices, corresponding to an infinity of analytic intcrpolants 
of the norming constants, see below) and thus arrive at a nontrivial jump across 
the two loops, encircling the segments [0,i^4] and [— iA, 0] respectively. This is a 
trivial deformation, valid for any h (not necessarily small). The discrete nature of 
the spectrum of L is mirrored in the discrete nature of the jump matrices: they 
involve a logarithmic integral with respect to a discrete measure. We sometimes 
refer to this as a discrete Riemann-Hilbert problem. 

THEOREM. Let dp = (p°(r)) + (p°)*(r]*))dr], where p° = i is the asymptotic 
density of eigenvalues supported on the linear segment [0, iA]. Set X(X) = Tr(X-iA). 

Letting M + and M_ denote the limits of M on £ = C U C* from left and right 
respectively, we define the Riemann-Hilbert factorization problem 

M+(A) =M_(X)J(X), 

where 

J(A) = v(X), X G C, 
= <7 2 «(A*)V2,AeC* ) 
limx^ocM(X) = I, 

and 
(8) 

u(A)= ^ J -iexp(\jlog(X-T 1 )d i i(r 1 ))exp(-\{2iXx + 2iX 2 t-X(X)))^ 

Then the solution of the initial value problem for the focusing NLS equation is 
given by tp(x,t) = 2i limx^oo(XMi 2 ). 

Note that in the statement of the theorem the measure in the logarithmic inte- 
gral is now continuous. We have effectively substituted a discrete set of eigenvalues 
by its continuous limiting density. This is only valid as h and the rigorous 
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justification of the discrete-to-continuum passage is far from trivial, especially near 
the point where the loops C, C* hit the eigenvalue spike. 

The analysis in [KMM] makes use of all the ideas described in the previous 
sections (factorization, lenses, the weighted Cauchy operator, an auxiliary scalar 
problem), but it also takes care of the fact that while the loops can be deformed 
anywhere away from the poles as long as h is not small, they have to be eventually 
located at a very specific position in order to asymptotically simplify the RH prob- 
lem, as h — > 0. Appropriately, the definition of a g-function has to be generalized. 
Not only will it introduce the division of the loop into arcs, called "bands" and 
" gaps" , but it must implicitly select a contour. Rather than giving the complicated 
set of equations and inequalities defining the g-function, we will rather focus on 
the associated variational problem; it is not a maximization problem but rather a 
maximin problem. Here's the setting. 

Let H = {z : Imz > 0} be the complex upper-half plane and H = {z : Imz > 
0} U {oo} be the closure of H. Let also K = {z : Imz > 0} \ {z : Rez = 0, < 
Imz < A}. In the closure of this space, K, we consider the points ix + and ix-, 
where < x < A as distinct. In other words, we cut a slit in the upper half-plane 
along the segment (0,iA) and distinguish between the two sides of the slit. The 
point infinity belongs to K, but not K. Define G(z; ry) to be the Green's function 
for the upper half-plane 

G{z , v) = log \^l 

and let dfj,°(j]) be the nonnegative measure — idr\ on the segment [0,i^4] oriented 
from to iA. The star denotes complex conjugation. Let the "external field" <f) be 
defined by 



G(z; vW (ri) - Re{n{iA - z) + 2i(zx + z 2 t)), 



where, without loss of generality x > 0. 
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Let M be the set of all positive Borel measures on K, such that both the free 
energy 

E(fi) = J J G(x,y)dn{x)dn{y), /i£M 
and J (fid/j are finite. Also, let 

V»{z) = J G{z,x)dn{x), (iel 

be the Green's potential of the measure [i. The weighted energy of the field <f) is 

E+in) = E(p) + 2 [ <f>dn, 



for any 

Now, given any curve F in K, the equilibrium measure A F supported in F is 
defined by 

E<f,(\ F ) = min li€M (F)E^(fj,), 

where M(F) is the set of measures in M which are supported in F, provided such 
a measure exists. 

It turns out that the finite gap ansatz is equivalent to the existence of a so- 
called S-curve joining the points 0+ and 0_ and lying entirely in K. By S-curve we 
mean an oriented curve F such that the equilibrium measure A F exists, its support 
consists of a finite union of analytic arcs and at any interior point of suppfi the so 
called S-property is satisfied 4 

d -{<t> + v xF ) = ^-^ + v^ F ), 



dn + dn. 
The appropriate variational problem is: seek a "continuum" 5 C such that 

(9) E^iX ) = max FeV £^,(A F ) = max FeV min^M^E^fi), 

where F is the set of continua lying in K. The existence of a nice S-curve follows from 
the existence of a continuum C maximizing the equilibrium measure, in particular 
the associated Euler-Lagrange equations and inequalities. 



^ ~r~ — , ~r~ — arc the normal outward derivatives on each side respectively 
5 a compact connected set containing 0+,0— 
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Problem (9) is the non-self-adjoint analogue of the Lax-Levermore problem and 
a nonlinear analogue of (3). 



4. JUSTIFICATION: EXISTENCE OF THE STEEPEST DESCENT PATH 

EXISTENCE THEOREM [KR]. For the external field <j>, there exists a contin- 
uum F 6 F such that the equilibrium measure A F exists and 

E <j> [F}(= E (t> (X F j) = max Fer min lie M(F)E ( t > {n). 

REGULARITY THEOREM [KR]. The continuum F is in fact an S-curve, so 
long as it does not touch the spike [0, iA] at more than a finite number of points. 

If F touches the spike [0, iA] at more than a finite number of points, a concep- 
tual revision is required. We briefly discuss this issue in the next section. 

Here are the main ideas of the proofs. Let po be the distance between compact 
sets E, F in K defined as 

Po{E,F) = max zeE min (eF p (z,Q. 

Introduce the Hausdorff metric on the set I(K) of closed non-empty subsets of K: 
Pk{A,B) = sup(p (A,B),p {B,A)). 

Compactness of F is the necessary first ingredient to prove existence of a max- 
imizing contour. The second ingredient is semicontinuity of the energy functional 
that takes a given continuum F to the equilibrium energy on this continuum: 

E : F —> E^[F] = E^(X F ) = itf p (E(p) + 2 J ^>dp). 

For regularity, the crucial step is 

THEOREM [KR] . Let F be the maximizing continuum of and X F be the equi- 
librium measure. Let x,t be such that F does not touch the spike [0, iA] at more 
than a finite number of points. Let p be the extension of \ F to the lower complex 
plane via p(z*) = —p(z). Then, if V is the logarithmic potential of p,, 

Re{ ( d J^± + V \z)f = Re(V'(z)f 2Re [ V ' {z) ~ V ' {u) dp{u) 
J u- z J z - u 

+Re[\ f 2{u + z)V\u) dp(u)}. 
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PROOF: By taking variations with respect to the equilibrium measure. 

It is now easy to see that the support of the equilibrium measure of the maxi- 
mizing continuum is characterized by 



The S-property follows easily and this proves the Regularity Theorem. 

It is worth mentioning here the recent work of Tovbis, Vcnakidcs and Zhou 
[TVZ] , which examines the initial value problem for the focusing NLS (in the semi- 
classical limit) under two different classes of initial data. Under one of these classes, 
no eigenvalues exist, hence the eigenvalue spike is missing. In such a case our argu- 
ment above would prove a regularity theorem without the extra assumption on the 
maximizing contour (that it does not touch the spike [0, iA] at more than a finite 
number of points). 



5. CROSSING THE EIGENVALUES BARRIER AND THE QUESTION OF 
SECONDARY CAUSTICS. 



If F touches the spike [0, iA] at more than a finite number of points, regularity 
cannot be proved as above because variations cannot be taken. In [KR] we have 
included a rough idea on how to extend the above proof. A more detailed argu- 
ment is forthcoming [K07] . Since a complete proof is not published yet we simply 
summarize our general plan. 

One wishes to somehow allow the contour F go through the spike [0,i^4]. One 
problem arising is that (the complexification of) the external field is not analytic 
across the segment [— iA, iA]. What is true, however, is that V is analytic in a 
Riemann surface consisting of infinitely many sheets, cut along the line segment 
[— iA, iA]. So, the appropriate underlying space for the (doubled up) variational 
problem should now be a non-compact Riemann surface, say L. Now, compactness 
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is the crucial element in the proof of a maximizing continuum. But we can indeed 
compactify the Riemann surface L by mapping it to a subset of the complex plane 
and compactifying the complex plane. The other problem, of course, is whether 
the amended variational problem (with the modified field defined on the Riemann 
surface and with the possibility of F not enclosing all the original eigenvalues) is 
still appropriate for the semiclassical NLS. The argument goes roughly as follows: 

(i) Proof of the existence of an S-curve F in L along the lines of [KR] . 

(ii) Deformation of the original discrete Riemann- Hilbert problem to the set F 
consisting of the projection of F to the complex plane. At first sight, it is clear that 
F may not encircle the spike [0, zA] . It is however possible to append S-loops (not 
necessarily with respect to the same branch of the external field) and end up with 
a sum of S-loops, such that the amended F does encircle the spike [0, iA\. To see 
this, suppose there is an open interval, say (ia,iai), which lies in the exterior of 
F, while ia, iot\ e F. Let us assume, that F crosses [0, iA] along bands at ia, ict\ 
(if not the situation is similar and simpler) ; call these bands 5*, Si . Let (3~ , [3 + be 
points (considered in C) lying on S to the left and right of ia respectively, and at a 
small distance from ia. Similarly, let /Sf,/?^ be points lying on Si to the left and 
right of ia\ respectively, and at a small distance from ia\. We will show that there 
exists a " gap" region including the preimages of fi~ , /3f lying in the Nth sheet for 
—N large enough, and similarly there exists a "gap" region including the preimages 
of j3 + , Pi lying in the Mth sheet for M large enough, both being regions for which 
the gap inequalities hold a priori, irrespectively of the actual S-curve, depending 
only on the external field! 

Indeed, note that the quantity Re{4> a (z)) (which defines the variational in- 
equalities) is a priori bounded above by — (f>(z). For this, see (8.8) in Chapter 8 of 
[KMM]; there is actually a sign error: the right formula is 

Re(4°(z)) = -<t>(z) + J G(z, r])p a (ri)dr]. 

Next note that the difference of the values of the function Re((j)' J {z)) in consecutive 
sheets is 5Re((j) a ) = ±2irRez, and hence the difference of the values at points on 
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consecutive sheets whose image under the projection to the complex plane is irj + e, 
where 77 is real and e is a small (negative or positive) real, is S(Re<p cr ) = ±2ire. 
This means that on the left (respectively right) side of the imaginary semiaxis, the 
inequality Re(<p' 7 (z)) < will be eventually (depending on the sheet) be valid at 
any given small distance to it. 

Applying the theory of [KR] we join the preimages of and (under the 
projection pf L to C) lying in the Nth sheet and the preimages of /3 + and (3f lying 
in the Mthe sheet. Finally we connect the preimage of f} + lying in the Nth. sheet 
to the preimage of (3~ lying in the Mth sheet, and so on. Note that the bands are 
S-curves with respect to any branch of the external field. We thus end up with an 
S-loop whose projection is covering the "lacuna" ia,ia\. 

The original discrete Riemann-Hilbert problem can be trivially deformed to a 
discrete Riemann-Hilbert on the resulting (projection of the) union of S-loops. All 
this is possible even in the case where F self-intersects. 

(iii) Deform the discrete Riemann-Hilbert problem to the continuous one with 
the right band/gap structure (on F; according to the projection of the equilib- 
rium measure on F), which is then explicitly solvable via theta functions. Both 
the discrete-to-continuous approximation and the opening of the lenses needed for 
this deformation arc justified as in [KMM] (see also the article [LM] mentioned 
below for the delicate study of the Riemann-Hilbert problem near the points where 
F crosses the spike). The g- function is defined by the same Thouless-type for- 
mula with respect to the equilibrium measure (cf. section 2 (iii)). It satisfies the 
same conditions as in [KMM] (measure reality and variational inequality) on bands 
(where the branch of the field turns out to be irrelevant) and on gaps (where the 
inequalities are satisfied according to the branch of the external field). 

We have sketched a proof that the solution of the initial value problem for the 
focusing NLS equation admits a global (in x, t) finite genus representation, asymp- 
totically as h — > 0, at least in the case of the simplest data u(x,0) = Asechx. But 
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the arguments above also hold for a large class of " semiclassical soliton ensembles" 
initial data defined precisely in section VIII of [KR]. 

It is worth noting here the recent paper [LM] which also adresses the issue of 
the target contour hitting the eigenvalue barrier. ([LM] does not really mention 
a variational problem and prefers to consider directly the conditions (equations 
and inequalities) for the g-function.) This very interesting paper does not prove 
the existence of an appropriate target contour but instead contains a numerical 
and theoretical discussion of the issue of the target contour hitting the eigenvalue 
barrier. 

In [LM] the "band" part of the contour is defined not as the support of an 
equilibrium measure but instead it is considered as the trajectory of a quadratic 
differential (as in [KMM]). It is noted numerically that such a trajectory may hit 
the barrier [0, iA}. It is then proposed that the inequalities defining the "gap" part 
of the target contour be amended when it passes the barrier [0, iA] and the actual 
amendment is justified numerically and theoretically. As a conclusion it is claimed 
that the mechanism of the second "caustic" (a caustic appears when the topology 
of F changes as x, t vary) is different from the mechanism of the first caustic. 

We simply note here that the extra conditions suggested in [LM] appear nat- 
urally in the framework we have introduced above. The different gap inequalities 
in [LM] correspond to gap inequalities in different sheets, as viewed from our per- 
spective. We thus prefer to say that the mechanism for any "caustic" (caused by 
the change of the topology of F as x, t vary) is essentially the same, independently 
of whether the maximizing contour has crossed the spike [0, iA] or not. In any case 
the real issue here is not whether we have a first or second or higher order caustic, 
but whether the maximizing contour has hit the spike [0,iA]. It so happens that 
for the very specific data Asechx the second caustic appears after the maximizing 
contour has crossed the spike [0, iA] once, but in general this does not have to be 
the case. 6 

6 Also, it is worth recalling Chapter 6 of [KMM] where it is noticed that, depending on the 
choices of the parameters a and J defined there, the maximizing contour may hit the barrier [0, iA] 
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6. CONCLUSION 

In the asymptotic analysis of Riemann-Hilbert problems arising from integrable 
systems where the associated Lax operator is non-self-adjoint, the computation of 
non-trivial steepest descent contours is essential. The two main components of a 
rigorous proof of asymptotic formulae are: 

(i) Proof of the existence and regularity of such steepest descent contours. 

(ii) Given (i), a rigorous proof of the asymptotic validity of the deformation of 
the given Riemann-Hilbert problem to one with jumps across the steepest descent 
contour. 

In this review paper we have presented some methods and results, contained 
in [KR] and [KMM], achieving (i) and (ii) for some specific cases of the initial 
value problem for the focusing integrable nonlinear Schrodingcr equation in the 
semiclassical limit. We expect that these methods and results may be useful in the 
treatment of Riemann-Hilbert problems arising in the analysis of general complex or 
normal random matrices [WZ] . Although there have been simple cases of non-self- 
adjoint problems where the target contour can be explictly computed without any 
recourse to a variational problem (which of course is always there; see e.g. [K96], 
[KSVW], [TVZ]), we believe that global results can in general only be justified by 
proving existence and regularity for a solution of a maximin variational problem in 
two dimensions. 
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